Momentum dependence of the superconducting gap in Bai_a;Ka,Fe2As2 
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The precise momentum dependence of the superconducting gap in the iron-arsenide superconductor 
with Tc = 32 K (BKFA) was determined from angle-resolved photoemission spectroscopy (ARPES) 
via fitting the distribution of the quasiparticle density to a model. The model incorporates finite 
lifetime and experimental resolution effects, as well as accounts for peculiarities of BKFA electronic 
structure. We have found that the value of the superconducting gap is practically the same for the 
inner F-barrel, X-pocket, and "blade" -pocket, and equals 9meV, while the gap on the outer F-barrel 
is estimated to be less than 4meV, resulting in 2A/A:bT'c = 6.8 for the large gap, and 2A/A;bTc < 3 
for the small gap. A large (77 ± 3%) non-superconducting component in the photoemission signal 
is observed below Tc. Details of gap extraction from ARPES data are discussed in Appendix. 



PACS numbers: 74.25.Jb 74.70.-b 79.60.-i 

I. INTRODUCTION 

Recently a new class of high-temperature superconduc- 
tors, iron-based pnictides, attracted much attention due 
to a rapid increase of the critical temperature, Tc, up to 
56 K.H These novel materials, still remaining terra incog- 
nita for theoreticians and experimentalists, require vast 
efforts from both sides to achieve a progress in the un- 
derstanding of their nature. One of the most important 
contributions that experimentalists can make to the de- 
velopment of a theory of any class of superconductors, 
is revealing the magnitude and symmetry of the super- 
conducting gap. Knowledge of the precise momentum 
dependence of the superconducting gap can provide de- 
sirable information about the pairing mechanism that 
underlies superconductivity in these compounds. Up 
to now there are a number of papers, providing differ- 
ent estimates of the superco nducting gap in iron-based 
superconductors,^ ^ l ^ l ^ l ^ l ^ l ^ l ^ l^ as well as different conclu- 
sions about the strength of coupling and applicability of 
BCS (Bardeen-Cooper-Schrieffer) theory to these com- 
pounds. Here we present an angle-resolved photoemis- 
sion spectroscopy (ARPES) study of the superconduct- 
ing gap in single crystals of Bai_a^Ka^Fe2As2, Tc=32K 
(BKFA) . The superconducting gap is extracted from pho- 
toemission data via a fit to a model, that accounts for 
finite self-energy, temperature, experimental resolution, 
as well as nonlinearity of the band dispersion, where it is 
necessary. 



II. RESULTS 

According to our recent study, the Fermi surface 
(FS) of BKFA, as seen in ARPES, consists of four differ- 
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Fig. 1 (color online), (a) Distribution of the photoemission 
intensity at the Fermi level (FL) with superimposed Fermi 
surface (FS) contours (white lines), (b) Momentum-energy 
cut through the F-point [cutl in panel (a)] taken at 10 K. 
(c) Same cut, taken at 45 K. (d), (e) MDC, taken at the FL, 
and symmetrized EDC from cuts (b) and (c) respectively. 
Maxima of the symmetrized EDC are marked by dots, (f) kp 
EDC referring to the inner F-barrel, recorded at 10 and 45 K. 
(g) Near-ZcF EDC emphasizes onset of the superconductivity 
even better, (h) Energy dependence of the inner F-barrel 
intensity, extracted from the fit of MDC. (i) The same for the 
outer F-barrel. 
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Fig. 2 (color online). Temperature and momentum depen- 
dence of the superconducting gap on the inner F-barrel. Re- 
liability of the results, (a) Evolution of the integrated EDC 
from cutl [see Fig. 1(a)] with temperature and fits to formula 
(1). EDC are shifted along left axis for clarity. Panel (b) 
shows the extracted magnitude of the gap, plotted against 
temperature. Inset to (b) shows temperature dependence 
of resistivity (with and without magnetic field), confirming 
high quality of the crystals and emphasizing equality of bulk 
and surface Tc = 32 K. (c) Integrated EDC from cut2 [see 
Fig. 1(a)], measured at different excitation energy at 11 K, 
and corresponding fits to formula (1) reveal reproducibility 
of the data and robustness of the fitting procedure. EDC 
are shifted along left axis. The inset in (c) shows a single 
EDC recorded with hu^SO eV, demonstrating high resolution 
at high excitation energies. 



ent sheets: outer F-barrel, inner F-barrel, X-pocket, and 
"blade" -pockets along the XF line^^ [see Fig. 1(a)]. X- 
pocket is electron-like, while all other FS sheets are hole- 
like. Fig. 1(b) and (c) show the same energy- momentum 
cut [cutl in Fig. 1(a)] through the distribution of the 
photoemission intensity at 10 and 45 K respectively. A 
backfolding dispersion of the inner F-barrel develops with 
cooling below Tc [Fig. 1(b), (c)], that points to the open- 
ing of the superconducting gap. To investigate the be- 
havior of the quasiparticle density near the Fermi level 
(FL) in detail, we plot symmetrized energy distribution 
curves (EDC) measured at 10 and 45 K in panels (d) and 
(e) respectively. The distance between the two peaks 
in the symmetrized EDC approximately equals to the 
doubled value of the energy gap, 2 A. As follows from 
Fig. 1(d) and (e), peaks in EDC, which correspond to 
the FL crossing of inner F-barrel, split into two below 
Tc, indicating the opening of a gap of the order of 9 meV, 
while peaks in the EDC, which refer to the outer F-barrel 
do not split upon cooling, indicating zero (or small in 
comparison with the peak width) magnitude of the gap 
on this part of the FS. Fig. 1(f), (g) show the energy 



dependence of the intensity [area under the momentum 
distribution curve (MDC)], which comes from inner and 
outer F-barrels, as extracted from the fit of MDC to four 
Lorentzians. A pile- up peak clearly develops on the curve 
that corresponds to the inner F-barrel, while no such fea- 
ture is observed for the outer F-barrel. The resolution- 
broadened 10 K Fermi cut-off is plotted in the panel (g) to 
show that the difference between 45 K- and 10 K-curves 
mainly comes from temperature smearing of the Fermi 
edge. The mentioned arguments allow us to conclude 
that the inner F-barrel bears a gap of the order of 9 meV, 
and the gap on the outer one is much smaller. 

Though straightforward and unpretentious, "sym- 
metrization" is a rough method for the gap extraction 
from photoemission data, therefore below we improve 
the assessment of the gap magnitude with a robust fit- 
ting procedure, where the value of the superconducting 
gap is extracted from the fit of EDC, integrated in a 
finite momentum window. In this case the integration 
is performed over a very small, compared to the Bril- 
louin zone size, region, which does not imply reduction 
to momentum-integrated data, and is used only in or- 
der to collect the whole available photoemission signal, 
referring to the particular FL crossing of a single band. 
The integrated EDC (lEDC) is fitted to the specially 
derived formula (see Appendix I), which coincides with 
Dynes function^^ multiplied by the Fermi function, and 
convolved with the response function: 
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where E = \/ {co -\- iYj"Y — A^, uj is the binding en- 
ergy with reversed sign, T is the temperature, T^" is the 
imaginary part of the self-energy, Ak is the momentum- 
dependent superconducting gap, and 5E is the exper- 
imental resolution. A similar method of gap extrac- 
tion is widely used in angle-integrated photoemission 
spectroscopy.^^ Fig. 2(a) shows lEDC, that refer to the 
inner F-barrel, from cutl [see Fig. 1(a)], measured with 
50 eV photon energy at different temperatures, as well as 
their fits to formula (1). The temperature dependence of 
the extracted gap, shown in Fig. 2(b), illustrates that a 
superconducting gap develops upon cooling through Tc, 
and reaches the value of 9.1 ± O.TmeV at low tempera- 
tures. Fig. 2(c) represents lEDC from cut2 [see Fig. 1(a)] 
recorded at IIK with different incident photon energies, 
hv. The data exhibit good reproducibility, and the values 
of the gap extracted for different hv show only a small 
scattering within the error bars — fit results in 9.4, 9.5, 
and 10.2 meV for hv = 40, 50, and 70 eV respectively. In 
order to emphasize the quality of our data recorded at 
high excitation energies, we show a single EDC recorded 
with hu=SOeV as an inset to Fig. 2(b). Thus, we can 
conclude that the momentum anisotropy of the super- 
conducting gap on the inner F-barrel is absent within 
1.5 meV. The outer F-barrel is much less intense than the 
neighboring inner F-barrel, which complicates the analy- 
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sis. With the same fitting procedure we estimate the gap 
on the outer F-barrel to be not more than 4 meV. 

Now we turn to the most interesting and problem- 
atic region of the BKFA Fermi surface, that was not 
completely resolved in previous studies of iron-arsenic 
superconductors — a propeller-like structure centered at 
the X-point [see Fig.l (a)]. Fig. 3(a) and (b) show the 
same energy-momentum cut through the X-point [cut3 
in Fig. 1(a)] above (36 K) and below (UK) respec- 
tively. Note that the intensity of the blades is largely 
suppressed for this photon energy and light polarization 
[see Fig. 1(c) in Ref. 11 . The difficulties with this region 
in momentum space are related to the presence of the van 
Hove singularity close to the FL, which brings the peak 
in the density of states already above [see Fig. 3(c)]. 
If both bottom and top of the band are far enough from 
the FL, then one can treat the dispersion of the band as 
linear without significant accuracy loss, so that formula 
(1) works well. For the case of band depth comparable to 
the magnitude of the superconducting gap, formula (1) 
has to be modified in order to account for the nonlinear- 
ity of the normal-state band dispersion. If one assumes 
that the band possesses electron-like parabolic dispersion 
with the bottom of the band located uo = — below 
the FL, then formula (1) transforms to 
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As it is easy to see, formula (2) reduces to (1) when sq 
becomes much larger than cj, E'', 5E and Ak. For a 
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Fig. 3 (color online). Superconducting gap on the X-pocket. 
(a) Energy-momentum cut through the X-point [cut3 on 
Fig. 1(a)] taken at 36 K. (b) Same cut taken at 11 K. (c) 
Evolution of the lEDC with temperature and fits to formula 
(2), and symmetrized /cf-EDC. (d) Shows temperature depen- 
dence of the gap. (e) Comparison of lEDC referring to the 
M-pocket and to the blades reveals virtually the same values 
of the superconducting gap. 




Fig. 4 (color online). Momentum dependence of the super- 
conducting gap in Bai-xKxFe2As2 (Tc = 32 K) is shown as a 
three-dimensional plot with underlying FS intensity map for 
orientation. Green line denotes the boundary of the Brillouin 
zone. 



detailed derivation, see Appendix I. 

The depth of the X-pocket was determined from the 
normal-state data using two different methods — a fit of 
momentum distribution curves taken at different binding 
energies with two Lorentzians, and a fit of the lEDC to 
formula (2) with A = and Sq as a free parameter. Band 
depths determined from both methods agree well — the 
first method results in sq = 20meV, while the second 
one yields Eq = 20.5 meV. Fig. 3(c) shows lEDC from 
cut3 [see Fig. 1(a)], referring to the X-pocket, measured 
at 11 and 36 K, as well as corresponding fits to formula 
(2). One may note in Fig. 3(c) the leading edge below 
the Fermi level for high temperature data, as well as two 
separate peaks in symmetrized EDO above Tc. These 
signatures of the gap are not relevant here, as discussed 
in Appendix II. Temperature dependence of the gap, ex- 
tracted from fitting the data to the formula (2), is shown 
in the Fig. 3(d). At low temperatures the gap on the 
X-pocket reaches 9.3 meV. From available experimental 
data we estimate the gap magnitude on the blades to be 
also 9meV [see e.g. Fig. 3(e)]. The results concerning 
momentum dependence of the superconducting gap are 
graphically summarized in Fig. 4. The gap is isotropic 
within the error bars, though, along with similarities to 
ReflU we see evidence for small anisotropy on the inner 
F-barrel — the gap may be slightly larger along FX (Bril- 
louin zone diagonal) than along FM (the difference is less 
than 10%). 

Presented analysis of the data via fitting of lEDC al- 
lows us to conclude that the low-temperature spectra 
have superconducting and non-superconducting compo- 
nents [see Fig. 5] . Only about 23 ± 3% of the intensity, 
coming from the inner F-barrel at 10 K, refer to the super- 
conducting part of the spectrum.^^ The presence of the 
two different components in the measured signal can be 
explained by a phase-separated coexistence of supercon- 
ducting and normal states, which was already observed 
in these^ as well as in other similar samples. 



3 



a) 



3 



±t 4 
CO 



C/) 

c 




b) 

3 100 

■D 



C/) 
C 





( 


O Experimental data yi 




— Fit % 




Non-superconductiriQ part 




I \ W 


Superconducting part \ 




1 1 1 1 
-40 -30 -20 -10 


1 1 

10 


CO (meV) 




Second derivatives 




Data 




Fit 




— SC part 




— Non-SC part 













c) 



lii 




-40 -30 -20 -10 10 

CO (meV) 

Fig. 5 (color online). Superconducting and non- 
superconducting constituents of the spectrum, (a) Energy 
distribution of the intensity, corresponding to superconduct- 
ing and non-superconducting parts of the spectrum, (b) Sec- 
ond derivatives of the data and fit. Structure of the sec- 
ond derivative confirms presence of superconducting and non- 
superconducting components, (c) Sketch, illustrating pres- 
ence of two different components in the same spectrum. 



III. EXPERIMENTAL DETAILS 

Single crystals of BKFA were grown using Sn as flux in 
a zirconia crucible. The crucible was sealed in a quartz 
ampoule filled with Ar and loaded into a box furnace. 
A cooling rate of 3C/h was applied from the maximum 
temperature of 850 down to 550 C for the growth. The 
growth details are described in Ref. 10, The crystals were 
cleaved in situ and measured with Scienta SES R4000 
analyzer at the base pressure of 5 • 10~^^ mBar. ARPES 
experiments were performed using the "1^ ARPES" end 
station at BESSY. 



Ref. num. HE] H [5] 

Tc 53K 37K 35K 53K 

Inner T-barrel 20 12.5 12 15 

Outer r-barrel — 5.5 8 — 

X-pocket — (12.5) (10) — 

Blades — — (11) — — 

Gap anisotropy — <3 2 <5<3 



m 

37 K 

12 

6 

(11) 



This paper 
32 K 
9.2 ± 1 

< 4 
9±2 
9±3 

< 1.5 



Table I: Momentum dependence of the superconducting gap in 
iron-arsenic superconductors, as revealed by ARPES studies, 
sorted by the time of appearance on the arXiv.org. Values 
of the gap and estimates of the gap anisotropy on the inner 
F-barrel are given in millielectron- volts. 



Ref. num. [2] [3] [4] [5] [6] [H 

Large gap 9 8.1 8.2 6.8 7.5 3.7 

SmaUgap — 3.6 5.5 — 3.9 — 



[8] |9] This paper 
9.6 4 6.8 

3.4 — <3 



Table IL Coupling strength, 2A//cbTc, in iron-arsenic su- 
perconductors, as revealed by different experimental tech- 
niques — compare to the BSC universal value 3.53. Most of 
the available studies reveal two superconducting gaps of differ- 
ent magnitudes, which are represented in the table as "large" 
and "small". Refs l2l3l4l5l6l are ARPES studies, Refs,7 8 are 
Andreev spectroscopy studies, Ref. 9 is a specific heat study. 



detection of the true values for the MDC width with 
an accuracy much better than momentum broadening. 
The lEDC-fitting procedure, applied to ARPES spectra 
of BKFA, yielded the following results: (i) the gap on 
the inner F-barrel along FM equals 9.1 ± 0.7 meV and 
along FX 9.7 ± 1 meV; (ii) the gap on the outer F-barrel 
is less than 4meV; (iii) the gap on the X-pocket equals 
9.3 ± 2meV; (iv) the gap on the blades is estimated to 
9 meV; (v) at 10 K the imaginary part of the self-energy, 
E'', in the vicinity of the FL was found to be equal 
to 1-2 meV. Comparison with other ARPES studies of 
the superconducting gap in iron-based superconductors 
is shown in Table L We evaluate the coupling strength 
as 2A//cbT'c=6.8 for the inner F-barrel, X-pocket, and 
blades, while for the outer F-barrel 2A//cB^ci3. A com- 
parison to other experiments is shown in Table IL 

Finally, the observation of drastically different super- 
conducting gaps on the inner and outer F-barrels is inline 
with theoretically suggested magnetic downfoldin^^^ and 
with a hidden (tt, 7r)-order observed experiment ally .H' 
Otherwise it would be hard to expect so different gaps on 
closely located and very similar bands, formed by slightly 
different combinations of the same atomic orbit als. 



IV. CONCLUSIONS 
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In conclusion, we have proposed a precise procedure 
for extracting the momentum dependence of the super- 
conducting gap from ARPES spectra. The developed 
method of data treatment allows to measure energy gaps 
with an accuracy much higher than experimental reso- 
lution, similarly to the Voigt-fit procedure,^ enabling 
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VI. APPENDIX I: FITTING FORMULA 



A. Derivation 

Below we adduce the detailed derivation of the formulae 
(1) and (2). We also show that for the simple case of neg- 
ligible curvature of the band dispersion and momentum 
independent gap, formula (1) coincides with the Dynes 
formula. 

A very general model for the measured ARPES signal 



By definition, the integrated EDC is 

lEDC(cj) = / I{k,uj)dk. 



Rk 



(3) 



(4) 



As soon as we anyway integrate our data over /c, momen- 
tum resolution does not affect lEDC,^^ which is already 
an advantage of this method. Substituting ([3| into Q, 
we get 



IEDC(a;) 



Ru^ (5) 



For the spectral function A(k^uo) in the superconduct- 
ing state, we use the following well accepted modeP^ 



A(k, uj) = 2^[ul5{uj - Ek) + vl5{uj + Ek)], (6) 



where 



/ 



A{k^uoQ)dk = 



1^0 1 



^sign(cjo) 



kfci I 



• (10) 



(In this formula ^ implicitly depend on cjq-) 

For the case of the linear band dispersion the derivative 
is constant, = const, and we arrive at 



/ 



A{k,ujQ)dk 



I wo I 



(11) 



This formula coincides with the Dynes function, al- 
though the premises for the latter are somewhat different, 
requiring the assumption of the momentum independent 
gap. Important difference in definition of our lEDC and 
well known Dynes function is that the former is a trace 
integral along one direction [see Fig. 6 (b)], while the lat- 
ter is a double integral over the whole momentum space 
[see Fig. 6 (c)]: 



Dynes (cj ) = j j A{\<i^uj) dk^dky. 



(12) 



Substituting here the aforementioned model for the spec- 
tral function (IgI), we go from a double integral to the 



Ek 



Ek 

Ek = \Kl + ^^- (7) 



Substituting ([6| under the integral in ([5| and omitting 
unnecessary constant factors, we get 



A{k^ ix)Q)dk 



1 

2 I wo 



Uf.5{ujQ-Ek) + v^5{ojQ + Ek] 
1 



1 ( J^ZR 

2 I wo 



dEk 



dfc 



k—ki 



k=k2 



dk 



(8) 




b) 



c) 





where ki^2 are solutions of (_ki,2 — ±\/^o ~ Below 
we will denote derivative by a prime: 
Expanding the derivative 



— 77/ 

dk \k=ki— ki' 



\El\ 



(9) 



we get 



Fig. 6 (color online), a) Spectral function in the supercon- 
ducting state for the case of linear normal-state dispersion 
^k- According to the formula (|6|, the spectral weight above 
the Fermi level is governed by Uf., and by below, b) Inte- 
gration in our case is performed along one energy-momentum 
cut (grey stroke), which intersects Fermi surface at only one 
point, c) In the momentum-integrated techniques integration 
is naturally performed over the whole momentum space. 
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integration along the contour 



Dynes (cc^o) = j) '^^l^^kl d/c 

^klV^kp^dA:. (13) 



When the depth of the band is much larger than the value 
of the superconducting gap, i.e. when we can neglect 
the nonlinearity of the dispersion (which is an important 
condition for the Dynes formula to hold!), this expression 
reduces to the integral over the Fermi surface: 



dk. 



(14) 



Dynes (cjo) = y |V£^k 

k:<ek=0 

Here we expand V^k similarly to formula (9), and get 

l^ol i_>. 1-1 



Dynes (cjo) 



k:^k=0 



|Va| dk. (15) 



As soon as coq and A (in this case) do not depend on k, 
one can pull them out from under the integral: 



Dynes (cjo) 



1^0 1 



^ |Va| 'dk. 



(16) 



k:a=0 



The integrand does not depend on cjq, therefore the whole 
integral is an unnecessary for our purposes constant fac- 
tor, which can be omitted, and we arrive at the same 
result as (11): 



Dynes (cjo) = 



1^0 1 



(17) 



B. Finite lifetime 



Up to now we have the result [formula (10)], obtained 
under the assumption of infinitely large lifetime, or, in 
other words, for very sharp bands. In such a case in or- 
der to get formula that incorporates effects of the finite 
lifetime, the following recipe is often used: take the for- 
mula, derived for infinite lifetime, add to the argument 
the imaginary part, and take real part of the result. 



g{u;) Re g{u;^iJ:''). 



(18) 



Below we show that in our case this trick provides the 
exact result. 

In order to account for lifetime broadening rigorously, 
one has to substitute the delta function in ([6| for a 
Lorentzian: 
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Fig. 7 (color online). Intergation along the contour on the 
complex plane. According to Cauchy's residue theorem, the 
integral along the contour Vr^t^ equals to the residue in the 
pole of the integrand inside, cjq + iT^' . 



which results in the possibility to rewrite the expression 
for the spectral function in the following way: 



J^" (k,ijj) = A(k,ijj)^ 



^L^'\i^-E^), 



(19) 



where A(/c, uj) stands for non-broadened spectral function 
(|6|. As convolution over uj commutes with integration 
over k. 



J A^" (k, u)dk = J A{k, u;)dk 



(20) 



We already know the result for integration of the spectral 



function over momentum — formula (10), and now the 
only problem is to evaluate the convolution. We will do 
it for linear and quadratic band dispersions, i.e. input 
parameters to derive formulae (1) and (2). 

Let g{uj) = jA{k,uj)dk, then in order to evaluate the 
convolution in (20), we have to calculate the integral 

+ CO 

/ g{ij)L^" {iJo - u;)du;. Res f{z) = 1 

-DO ^ = ^0 

The function g{u;) is defined on the real axis. Once 
we know the analytic function g{z)^ 2: G C, such that 
Re{g{uj)) = g{uj) for G M, we can calculate the required 
integral with the help of Cauchy's residue theorem: 



—00 

J g{ijj)Lp {ujQ — uj)duj = 

-co 

— CO 

J g{u;)L^ {uoo - cj)d 



+ CO 



Re 



Re 



lim (b g{z)L (cjq — z)dz 

R^oo,r]^0 J 

I g possesses no poles inside F^^^^ I 
27rz • Res g[z)- 



z=uoTi^" 27r (cjo - ^)2 + S'^^ _ 

Regiooo^i^''), (21) 
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which coincides with formula (18), and imphes formulae 
(1) and (2) as corollaries. For definition of the integration 
contour Tji^^^ refer to Fig. 7. Explicit form of the function 
g{z) for linear band dispersion is 

where for z = re^^ we pick the following definition of the 
square root ^ = ri/2gi0/2^ ^ ^ [0,27r). 
For quadratic dispersion we get 



1 z /I 1 \ 1 / 1 1 



2 ^/z^ - A2 

where A:i,2 = V^o ± a/^^^^, = ^i/2gi0/2^ ^ ^ 
[-7r,7r). 

Defined in such way, ^1,2 (^) are analytic in C \ 
(—00, +00), i.e. all conditions for Cauchy's residue theo- 
rem are fulfilled. 

C. Formulae in real numbers 

For numerical calculations it is useful to rewrite for- 
mula (1) without the use of complex numbers: 



lEDC(cj) 



uj{a + c) + T."h 



\f2c\J~a- 



where a = cu^ - T."^ - Al, b = 2S''cj, and c = Va^TF. 
Similarly, formula (2) can be rewritten as 



1 



/(c^,T)--V^- 



lEDC(cj) = 

\u;[{a + c){ai + 71) + bf3i] + ^''[b{ai + 71) - A(a + c) 
2c7i + c^ai + 71 

|cj[(a + c)(a2 + 72) + ^/^2] + ^''[b{a2 + 72) - /52(a + c) 



2c72 + c^a2 + 72 



• sign(cc;) 



V^i + 71 V<^2 + 72 



^271 V272 
where 

a = 00^- Y."'^ -/Si, b = 2Y!'uj, c = Va^^b^, 

b 



<^i,2 = £0 T \/ -^^1 Pi, 2 = ± 



and 71,2 = \/«i,2 +/^i,2- 



VII. APPENDIX II: EXTRACTION OF THE 
GAP FROM THE MODELED DATA BY 
"SYMMETRIZATION", "LEADING EDGE", AND 
FITTING 





"Symmetri- 


Leading 


Fit 




zation" 


edge shift 







10 


9.1 


10 ±0.1 


4 


10 


8.6 


10 


8 


10 


8.0 


10 


12 


9.8 


7.8 


10 


16 


8.6 


7.5 


10 


20 





7.1 


10 ±0.1 



Table IIL Superconducting gap, as extracted from modeled 
data (Fig. 8) with different methods. All numbers are given in 
millielectron- volts. The correct value of the gap (implemented 
in simulation) equals lOmeV. 



metrization" is to be substituted by more rigorous and 
advanced ways of data treatment, such as those used in a 
very recent publications on photoemission spectroscopy 
of superconductors.^^ 

Below we model ARPES spectra with formulae ([3| and 
Q, and extract the gap with "symmetrization" , "lead- 
ing edge",^^ and proposed here fit of the lEDC. Results 
confirm that the fitting procedure is robust against mo- 
mentum resolution, properly accounts for energy resolu- 
tion and finite lifetime, provides correct values even in 
the case of the nonlinear band dispersion, and allows one 
to disentangle non-superconducting and superconducting 
parts of the spectrum. At the same time, "symmetriza- 
tion" and "leading edge" are not stable with respect to 
the effects of the experimental resolution, and further- 
more fail in the case of the shallow band and in the pres- 
ence of the non-superconducting component. 



A. Energy resolution 

First, we study the influence of the experimental en- 
ergy resolution on the determination of the gap from 
ARPES data with "symmetrization" , "leading edge", 
and fit to formula (1) from the Manuscript. The results 
of these studies are shown in Fig. 8 and summarized in 
Table III. Please note that not the resolution of the an- 
alyzer is important, but the resolution of the whole pho- 
toemission experiment. Also it is worthwhile to mention 
that effects of the lifetime broadening are in some respect 
similar to the effects of energy resolution, as they both 
lead to the broadening of the spectra. 

By the way, leading edge (the lowest binding energy at 
which the kp EDC reaches half of its maximum) alone 
is not a good measure of the gap (see corresponding 
columns in Figs. 8-10), while leading edge shift (shift of 
the leading edge with respect to the position in the nor- 
mal state) is a lot more relevant quantity. 



Momentum resolution 



The "symmetrization" is highly valued by some part of 
the ARPES community. We strongly believe that "sym- 



Next, we consider the influence of the experimen- 
tal momentum resolution on the determination of the 
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Fig. 8 (color online). Influence of the energy resolution on the determination of the gap via symmetrization, leading edge, and 
fit. First column: energy resolution for the corresponding row. Second column: simulated energy- momentum cut above Tc 
ij^" — 3meV, kT =3meV, A — OmeV). Third column: simulated energy-momentum cut below Tc ij^' — 3meV, kT =lmeV, 
A = lOmeV). Fourth column: determination of the gap with "symmetrization". Fifth column: determination of the gap with 
"leading edge". Sixth column: determination of the gap with fit to formula (1) and criterion as insets to some panels. 
"Symmetrization" and "leading edge" provide acceptable results for good resolution, and fail when the resolution becomes 
worse. The fitting procedure always provides the correct result. 
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"Symmetri- 


Leading 


Fit 




lO-^A"' 


zation" 


edge shift 










10.0 


9.1 


10 





50 


12.0 


8.2 


10 





100 


12.3 


8.0 


10 





200 


12.4 


7.5 


10 


10 


100 


16 


6.4 


10 


20 


200 


23 


4.9 


10 ±0.5 



Table IV: Superconducting gap, as extracted from modeled 
data (Fig. 9) with different methods. All numbers, except for 
momentum resolution, are given in millielectron- volts. The 
correct value of the gap (implemented in simulation) equals 
lOmeV. Parameters Ruj — lOmeV and Rk = 0.1 A corre- 
spond to the experimentally observed widths of the spectra. 



Rk, 


Input to 
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Fit 




L the model 


zation" 


edge 













-1.9 


+ 0.8 




10 


10.1 


6.9 


10 ±0.1 


45 





10.9 


0.5 


0±0.8 




10 


12.4 


7.9 


10 ±0.1 


60 





18 


3.5 


0±0.8 




10 


20 


8.3 


10 ±0.1 



Table V: Superconducting gap, as extracted from modeled 
data (Fig. 10) with different methods. All numbers, except 
for momentum resolution, are given in millielectron- volts. 



gap from ARPES data with "symmetrization" , "lead- 
ing edge", and the fitting to formula (1). The results 
of these studies are shown in Fig. 9 and summarized in 
Table IV. Note that not the resolution of the analyzer is 
important, but the resolution of the whole photoemission 
experiment. 

The width (full width at half maximum) of the nar- 
rowest EDC from Refs. 2-6, as well as from our studies is 
8-10 meV. The momentum resolution is about 0.1 A 



C. Nonlinearity of the band dispersion 

The case when the band depth is comparable to the 
value of the superconducting gap is quite complicated, 
and really requires special treatment. That is why for- 
mula (2) has been derived and used to fit the data. It is 
easy to mistake the van Hove singularity for the gap when 
using simplified methods of data analysis. Masking ef- 
fects of van Hove singularity is one of real examples where 
"symmetrization" and "leading edge" give wrong results 
(Fig. 10, and especially column three, bottom row). Nat- 
urally, such "gap" will not close at Tc. 

Here we have modeled the influence of the nonlinearity 
of the band dispersion in conjunction with experimen- 
tal momentum resolution on the determination of the 
gap from ARPES data with "symmetrization" , "leading 
edge", and fit to formula (2). The results of these studies 
are shown in Fig. 10 and summarized in Table V. 



D. Non-superconducting component 

According to recent /iSR (muon spin rotation) 
studies, superconducting fraction for optimally doped 
Bai_a^Ka^Fe2As2 samples (those used in Refs. 2-6) com- 
prises 50% of the sample volume (see Ref.|T7j), and for 
our slightly underdoped samples it comprises only 25% 
(see Ref. 16). Under these circumstances, "leading edge" 
is completely irrelevant to the gap value, while "sym- 
metrization" may provide some estimates of the value of 
the gap depending on other conditions (resolution, life- 
time etc.). Fitting in this case is indispensable, as it not 
only reveals precise values of the gap, but also allows 
to determine the fractions of the superconducting and 
non-superconducting signals (see Fig. 5). For our crys- 
tals these fractions determined from two completely dif- 
ferent methods — /iSR and the fit of the ARPES data — 
perfectly match each other. 



E. Renormalization 

Presence of the dispersion anomalies, "kinks" can af- 
fect position of the leading edge and peaks in the sym- 
metrized EDC, and can be mistaken for energy gap, sim- 
ilarly to the discussed above van Hove singularity. In 
the case of linear bare band dispersion,^ lEDC is not 
affected by self energy at all: 



+ CO 



1 



27r {co - ^'{uj) - v^kf + T."{uoy 



-dk 



1 



(24) 



F. Summary 

1. The proposed fitting procedure is rigorous and pre- 
cise method of gap extraction, which accounts for 
several important features of Bai_a;Ka^Fe2 As2 pho- 
toemission spectra: 

(a) nonlinearity of the band dispersion; 

(b) presence of large non-superconducting compo- 
nent; 

(c) experimental resolution. 

2. "Symmetrization" is not a universal way for the 
extraction of the gap from spectroscopic data, since 
it is highly sensitive to experimental resolution, and 
non-linearity of the band dispersion. For example, 
it 

(a) gives zero value for the gap while there is sub- 
stantial gap (Fig. 8, bottom row); 

(b) gives substantial value for the gap, while ac- 
tual gap is zero (Fig. 10, bottom row). 



9 



meV 10VA 



50 



100 g 



200 



10 100 



20 200 



above 7, 



below T, 



symmetrization 



leading edge 




k{Mk) 



(D (meV) 



CD (meV) 



(D (meV) 



Fig. 9 (color online). Influence of the momentum resolution on the determination of the gap via symmetrization, leading 
edge, and fit. First column: energy resolution for the corresponding row. Second column: momentum resolution for the 
corresponding row. Third column: simulated energy-momentum cut above Tc {T^" = 3meV, kT = 3meV, A = OmeV). Fourth 
column: simulated energy- momentum cut below Tc {T^" — 3meV, kT — 1 meV, A = lOmeV). Fifth column: determination of 
the gap with "symmetrization" . Sixth column: determination of the gap with "leading edge" . Seventh column: determination 
of the gap with fit to formula (1) and criterion as insets to some panels. "Symmetrization" and "leading edge" provide 
acceptable results for good resolution, and fail when the resolution becomes worse. The fitting procedure always provides the 
correct result. 



3. "Leading edge" alone is not a good measure 
of the gap (see bottom rows of the Figs. 8-10), 
while leading edge shift in absence of the non- 
superconducting component is a quite good, al- 



though still rough measure of the gap, and provides 
result with an accuracy better than 50% even under 
severe conditions (see Tables III-V). 



10 




Fig. 10 (color online). Influence of the small band depth on the determination of the gap via symmetrization, leading edge, and 
fit. First column: momentum resolution for the corresponding row. Second column: simulated energy- momentum cut above Tc 
{Y^" — 3meV, kT — 3meV, A — OmeV, Eq — 20meV). Third column: simulated energy-momentum cut below Tc {T^" — 3meV, 
kT — 1 meV, A = lOmeV, £q — 20meV). Fourth column: determination of the gap with "symmetrization". Fifth column: 
determination of the gap with "leading edge" . Sixth column: determination of the gap with fit to formula (2) and criterion 
as insets to some panels. First row: no resolution effects added. Second row: small resolution effects are added. Third row: 
moderate resolution effects are added (resolution effects are comparable, and may be even smaller to those in real data, which is 
easy to see comparing these energy-momentum cuts to directly measured). For simplicity, only momentum resolution is added. 
"Symmetrization" and "leading edge" provide acceptable results for good resolution, and fail when the resolution becomes 
worse. The fitting procedure always provides the correct result. 
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fi dx • 




which in our case reads 




12 



